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ABSTRACT 


This paper discusses several approaches to scattering from 
slightly rough, very rough, composite, and rough spherical surfaces. 

Average incoherent scattering cross sections for slightly rough 
surfaces are obtained using a perturbation technique for a perfect con- 
ductor, and for a surface material which is homogeneous (with variable 
permeability, p r , as well as permittivity, €r) * The results, correct 
to the first order, include polarization dependence, and particular 
solutions with curves are presented for the vertical, horizontal, and 
circular states. They are compared with measured data. The scatter- 
ing process for this class of surfaces is interpreted physically. 

Three optics techniques yielding results for very rough surface 
scattering cross sections are presented, and it is shown that they all 
give the same solution for Gaussian surfaces. Both homogeneous and 
perfectly conducting surface materials are treated, and polarization 
dependence is preserved. The restrictions under which the solutions 
are valid are enumerated. The physical interpretations of the scatter- 
ing mechanism for this class of surfaces are discussed and shed valuable 
insight on the process. Curves are shown for backscattering, employing 
two probability models . 

Surfaces consisting of both very rough and slightly rough 
structures together, termed composite surfaces, are analyzed, and an 
heuristic physical derivation of their scattering cross section is pre- 
sented. Curves for backscattering cross section are shown, and com- 
parison with measured results confirms the validity of this approach. 

Average backscattering cross sections for rough spherical sur- 
faces are presented. Both the coherent and incoherent cross sections 
are given for slightly and very rough spheres. These results have 
application to planetary surface radar scatter and to passive satellite 
communications . 

The large volume of published material on the subject is reviewed. 
The first two parts dealing with slightly rough and very rough surfaces, 
rather than retrace detailed mathematical derivations available elsewhere, 
attempt to accomplish five goals: (i) to compare the various analyses 
already available, (ii) to elucidate explicitly the approximations involved 
in each approach, (iii) to explain and interpret the physical process be- 
hind the mathematics producing the scattering, (iv) to present meaningful 
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closed form mathematical solutions and curves, many of which have not 
appeared before, for the scattering cross sections, and (v) to compare 
theory with measurements. The material in the last two sections, 
treating composite and rough spherical surfaces, is new as far as the 
Western literature is concerned, and is based directly on the results of 
the first two sections; physical interpretation of the scattering mechanism 
again is emphasized. 
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SCATTERING FROM SURFACES WITH DIFFERENT 
ROUGHNESS SCALES; ANALYSIS AND INTERPRETATION 

I. INTRODUCTION 

Serious attempts to analyze the scattering of electromagnetic 
waves by rough surfaces and interfaces did not begin until about 1950. 
Research on this subject was stimulated to a large extent by an ac- 
cumulation of measured data on radar scattering from terrain and the 
sea which was made during and immediately after World War II. 

Another stimulating factor was the significant expansion and dis- 
semination of knowledge in the area of applied statistics which took 
place in the decade before 1950. 

Since then, analysis of rough surface scattering has made use of 
one of two classes of surface model. The first, which we shall term 
the geometrical model, deals with surfaces made up of given determin- 
istic shapes, but arranged or distributed in a random fashion; distri- 
butions of hemispherical bosses on a plane is an example of such a 
model analyzed by Twer sky (1957). The second which we shall call 
the statistical model, treats the roughness height itself as a random 
variable from the outset. In the latter case, one must choose or specify 
the surface probability distributions or correlation coefficient. The 
statistical models seem more physically reasonable, since nature 
rarely composes a surface of given deterministic shapes. It is this 
second class of rough surface models which concerns us here. A 
review of the literature covering many of the rough surface scattering 
investigations is given in each respective section. More details, and 
results on other geometrical and statistical models can be found in 
Radar Cross Section Handbook (1968) . 

Statistically rough surfaces can be divided into three classes, 
distinguished by the relative roughness height compared to wavelength: 
slightly rough surfaces, intermediate, and very rough surfaces. The 
intermediate class has roughness features which lie in a region where 
neither high nor low frequency approximations apply. Only the slightly 
rough and very rough classes can be handled mathematically. In 
general, different techniques must be employed to solve these two 
classes. Much information can be obtained about the scattering 
behavior of rough surfaces from a study of these two classes, and 
it is not unreasonable to expect that extrapolation between them can 
give valuable insight into scatter from the intermediate class of rough 
surfaces. 
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Unfortunately, as with all new fields of investigation, errors and 
oversights have occurred. Soviet and Western developments in this 
area often paralleled each other; most investigators 1 articles and 
results went unnoticed by their foreign counterparts. Frequently, 
controversy in the literature over mistakes or minor details tended 
to obscure the positive and common points of agreement. Most 
serious of all, the emphasis on mathematical detail has in the 
majority of cases tended to obscure the simple physical concepts 
and processes producing the scattering. Also, many invetigators 
were not careful in stating the approximations and assumptions 
inherent in their analyses; this has often led to misinterpretation 
and misuse by others . 

We shall attempt here to emphasize the physical explanation 
of the scattering mechanism behind the mathematical formulations 
and solutions . Care is taken to outline the approximations and 
restrictions which apply to each model. Nonetheless results are 
presented in the most general form possible, including the effects 
of polarization and surface constitutive parameters . Several new 
results and explanations are offered for slightly rough surfaces, 
surfaces consisting of composite roughnesses of several scales, and 
rough spherical surfaces. In the case of very rough surfaces several 
existing analyses based upon different physical formulations are 
compared and shown to give identical results. 

To introduce the terminology we shall use, imagine that a 
perfectly smooth surface becomes slightly rough. The average 
specularly reflected field present for the smooth surface decreases 
somewhat as the roughness grows; in addition there is a smaller 
amount of power scattered away from the specular direction due to 
the presence of the roughness. The former we shall call the coherent 
component while the latter is termed the incoherent component of the 
scattered field. More precisely, if one surface of an ensemble is 
replaced by another statistically similar member, or if the surface 
is translating slowly with respect to transmitter and receiver, the 
% coherent component is distinguished by the fact that the average value 
of the complex field (existing at and near the specular direction) is 
non-zero. The incoherent field fluctuates, however; its phase angle 
is uniformly distributed; its average value is zero. On the other hand, 
the average scattered incoherent power (or incoherent field magnitude 
squared) is non- zero. As the roughness size increases, the coherent 
component decreases and the incoherent component increases. In the 
very rough surface limit, all of the scattered power is incoherent 
although it may be sharply beamed in the specular direction. This paper 
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is concerned mainly with the incoherent scattered power, except for 
the brief discussion below. 

The coherent field scattered from a slightly rough surface is 
computed in the same fashion as it is for a smooth surface. In con- 
trast to the incoherent component, it is highly dependent upon the 
shape of the illuminated area. It exhibits the familiar lobed structure 
of a flat plate, for example, while the incoherent component does not. 
Being concentrated around the specular direction, it can often be ignored 
in other directions. Experimentally, one can isolate the component by 
averaging first and then squaring the received voltage. The specular 
reflection coefficients for the slightly rough surface may be written in 
terms of the Fresnel coefficients for a smooth surface and a factor 
which accounts for roughness. 

This modification is discussed in Radar Cross Section Handbook 
(1968) and can be determined for a Gaussian surface height. The result- 
ing electric field reflection coefficients for that case are 


( 1 ) 


R 


R 11 (®i) = Rm (©i) e 


-2ko h. 2 c o s 2 9 j 


RfOi) = R x (6i) e 


where R M (0^) and R x (0^) are the usual Fresnel reflection coefficients for 
the smooth surface and are given in Eq. (35) . The quantity 9^ is the 
angle of incidence from normal, k Q = 2tt/\ is the wave number, and h 
is the rms roughness height. As can be seen from the above, the effect 
of roughness on the coherent scattered field is contained in the ex- 
ponential factor. It is this factor which quantitatively expresses the 
well-known Rayleigh criterion. Clearly, for k Q h < l/4, the coherent 
component will be not significantly different from that of a smooth surface, 
and we call the surface "slightly rough". On the other hand if k 0 h cos 9j> 1, 
the coherent component will be insignificant, and we call the surface 
"very rough". The coherent component disappears as the roughness 
height, h, increases. 

II. SLIGHTLY ROUGH PLANAR SURFACE 
A. Introduction 

In this section we will discuss an important class of random rough 
surfaces: those having a small scale of roughness whose rms roughness 
height, h, is much less than wavelength. In particular we shall obtain 
expressions for the incoherent average scattering cross sections per 
unit surface area, cr°, for various polarization combinations. As 
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mentioned in the introduction, this cross section accounts for most of 
the scattered energy everywhere except in the neighborhood of the 
specular direction. At and near this direction, the coherent component 
is strong and cannot be neglected. 

Slightly rough surface models have been employed in the past 
based upon the tangent plane, or Kirchhoff approximation . The first 
such analysis was that of Davies (1954). It has been discussed also by 
other workers, and a readable treatment is found in Beckmann and 
Spizzichino (1953 Section 5.3). It should be mentioned, however, that 
those models are often inapplicable to natural surfaces. The reason is 
that the tangent plane approximation restricts one to surfaces having 
ZiTp/X. > > 1, where p is the local radius of curvature of the surface 
at all points. This requirement, along with the defining requirement 
ZttIi/ X < < 1 generally describes a very smooth, undulating surface with 
hills much less than wavelength in height and hundreds of wavelengths 
apart; the maximum slopes of these hills are infinitesimal. What is 
more, the model fails in the low frequency limit, since wavelength 
eventually becomes larger than surface radii of curvature. Also, the 
model shows no polarization dependence for backscatter, as is typical 
of models based upon the tangent plane approximation. Experimental 
evidence shows that radar backscatter does depend upon polarization; 
this will be discussed again later. 

For the above reasons, we shall avoid use of the Kirchhoff 
approximation for the slightly rough surface model. Instead, we will 
resort to a much less restrictive model, but interestingly enough, a 
model originally formulated before that of Davies by Rice (1951). 

Rice's formulation was later developed by Peake (1959) who derived 
the scattering cross-sections. This model is based on a perturbation 
approach. Since the tangent plane assumption is not employed, all 
of the accompanying approximations such as neglect of multiple scat- 
tering and shadowing are avoided. The model is therefore valid for 
grazing incidence and scattering angles, and exhibits meaningful 
polarization dependence. Most important, it is valid in the low fre- 
quency limit and can apply to any of the class of very rough surfaces 
to be examined later when the frequency is low enough. 

Valenzuela (1967) recently extended the same technique to obtain 
expressions for the second order perturbation correction terms for 
backscatter from dielectric and perfectly conducting surface. Only the 
first order terms of the solution are obtained here . 

The following are the restrictions on the class of surfaces to which 
the model applies: 
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(a) k Q 4(x, y) < 1.0, i. e, , roughness height is small compared 
to wavelength, k; £>(x, y) is the height of the rough surface above the 
x-y plane, which is taken to be the mean surface plane. 

(b) 9£>/ 3x, 9£/9y < 1.0, i.e. , surface slopes are relatively 

small. 

(c) <3C,/3x 2 > = < 3^/3y 2 >, i.e., the roughness is isotropic. 

Here, the brackets < > indicate average over an ensemble of surfaces. 
This restriction is not essential to the solution, but is employed here 
for simplicity. 

(d) L >> i,\, i.e., the dimension, L, of the illuminated area 

is much greater than either the roughness correlation length, ft , or 
the wavelength. 

Since the tangent plane approximation is not employed, the mutual 
interaction effects of shadowing and multiple scattering are not over- 
looked. Hence the slightly rough surface solution includes these 
mechanisms. In fact, mutual interaction between neighboring hills 
is an important contributor when the distance between them is less 
than a wavelength. 

We shall obtain first results for the average incoherent scat- 
tering cross section per unit surface area for a surface of homogeneous 
material having relative permeability, p r , and permittivity, € r . In this 
respect our results differ from the analysis of Rice (1951), Peake (1959), 
and Valenzuela (1967) , who treat only dielectric surfaces. The 
quantities p r and e r may be real or complex, representing a lossless 
or lossy material. Their magnitudes may be greater, equal to, or 
less than unity, so that the important class of plasma media can be 
included. As a special case, the imaginary part of c r can approach 
infinity, yielding the proper results for a perfectly conducting surface. 
The notation of Rice (1951), Peake (1959), and Valenzuela (1967), will 


*Also implied in the definition of scattering cross section is the restriction 
that the observation point must be in the far field of the surface. In the 
absence of roughness, this requirement means that the distance to the 
observation point, R, must satisfy the criterion R > 2L» 2 /\. However, 
it can be shown (see Barrick, 1965) that when one is considering the 
incoherent scattered power from a roughened surface, this requirement 
reduces to R > 2f 2 A* where l is the roughness correlation length. 

Physically this means that one can be considerably closer to a rough 
surface than a smooth one, since ft < < L by (d) above. 
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be employed as much as possible here, so that readers wishing a more 
detailed derivation can refer to that work. The solution for a horizintally 
polarized incident plane wave will be obtained here. Similar results 
for vertical incident polarization will be then written from inspection 
by interchanging the roles of E and H, <= r and p r . From these results, 
the pertinent values for the circular polarization states and for linear 
polarization states of arbitrary orientation will be given. 

B. Analysis 

The time convention used here will be e’ la3t . The analyses of 
Rice, Peake, and Valenzuela are compatible with this if one replaces 
i there by -i, since they use the positive time convention. Let us first 
consider one particular rough surface from an ensemble of statistically 
similar rough surfaces. For convenience, assume the surface to be 
square, of length L, centered about the origin and with its mean surface 
coinciding with the x-y plane (Fig. 1). Later in the analysis, the 
surface length will be permitted to approach infinity; the square shape 
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of the surface is irrelevant to the results, but is a useful artifice in 
performing the derivation. Under these conditions, the surface height 
can be expanded in the following Fourier series: 

CO 

( 2 ) 

m, n=-oo 

, 2tt 

where a = — . 

L 

The real nature of the surface height requires that P(-m, -n) = 

P* (m, n) where the asterisk denotes complex conjugate. The trick in 
simplifying the analysis is to expand the scattered, or perturbed, field 
into a Fourier series with eigenfunctions e ^ a ( mx "^ n y) 0 Physically this 
means representing these fields as a superposition of plane waves in 
all directions. Remembering that the incident field is horizontally 
polarized in the y-direction as shown in Fig. 1, the total E -field in 
free space above the surface is then written as follows: 


(3a) Ej = Ej 

(3b) Ey = E 1 + E r + E®, 

(3c) E+ = E|, 


where E 1 is the incident field, and E r is the reflected field from a 
perfectly smooth homogeneous surface. The terms with the super- 
script, s, are then due solely to the presence of roughness, and may 
be considered the incoherent scattered or perturbed fields. These 
quantities are given as 


(4a) 


E 1 = E°e- aax+<rYZ ; E r = REje^^ 2 , 


(4b) 


E I = 



. ia(mx+ny ) ib(m, n) z 
Amn e c * 


m, n- 
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(4c) 


(4d) 



CO 


Bmn 

m, 


e ia(mx+ny) 


ib(m, n) z 

® a 



c 


mn 


e ia(mx+ny) 


m, n = 


e ib(m, n) z 


In these equations, cr = -ik Q = iu>J e 0 p 0 , a = sin 0j, Y = cos 0^, 
and E° is the electric field strength of the incident field. R is the 
Fresnel reflection coefficient for horizontal polarization, and is 
given by 


(4e) 


1 - 


R = 


t Y' 
p r o-Y 


1 + 


t Y' 


where T = 


p o-Y 


= ik o/VV and Y ' =y= 


r^r 


a 


The function b(m, n) is not independent of the wavenumbers in the 
x and y direction, and is dictated by the requirement that each term of 
the above series must satisfy the wave equation. This results in the 
following relationship: 


(5) 


b(m, n) 


= J k o - 


a (m 2 + n 2 ) 




2 . 2 , 2 . 2 , 
cr + a (m + n ) M 


The objective is to derive A mn , B mn , and C mn in terms of 
P(m, n) and all of the above parameters. This is done by a standard 
perturbation technique and the details are similar to those of Rice, 
Peake, and Valenzuela. Briefly, these three coefficients are obtained 
by expanding all quantities possible into series, each term of which 
involves a higher power of the smallness parameters, which are taken 
to be k 0 t,, 9t,/3x, and 9^/9y. By assumptions (a) and (b) , these 
quantities are small, and consequently only the first term of these 
series need be retained. The above three unknown coefficients (as 
well as three other unknown coefficients, G mn , H mn , I m n 
represent the scattered fields inside the surface material) are 
obtained from six equations, four of which come from the boundary 
conditions for the E and H fields at the interface z = and two of 
which arise from the divergence conditions on each side of the interface. 
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The three desired coefficients A mn , B mn , and C mn are then 
correct to the first order, i. e„ , the neglected terms are at least of 
order k 0 £,, 3£,/3x, or 3^/9y with respect to the result shown. The 
terms of the next higher order are derived by Rice (1951) and 
Valenzuela (1967) for perfectly conducting and dielectric surfaces, 
but they are not given here. Also it should be noted that each of 
these coefficients is directly proportional to P(m-v a n), the m-v, 
n-th coefficient in the Fourier expansion for the surface itself. 

(v = k Q sin 0p/a, and is taken to be an integer here, i. e. , the 
possible directions of incidence are quantized for the moment. 

However, in the limit as a = 2 tt/L — ► 0 with L co , v becomes 
large, and thus it is possible to approximate it by the nearest integer 
with negligible error, ) 

Up to now, no restrictions have been placed on the nature of 
the surface. The surface may in fact be deterministic or periodic. 

For example in the case of a one -dimensional sinusoidal surface of 
period L, the expansion for the surface, £,(x) , degenerates to the 
situation where P(j, k) = 0 for j ^ - 1 and k 0; this important 
case was solved long ago by Lord Rayleigh (Dover, 1945). To gain 
insight into the scattering process for the slightly rough surface, it 
is instructive to point out a fact about this surface. The vanishing 
of all but two of the coefficients P(m-v, n) means that the coefficients 
for the scattered field expansions of Eqs. (3b), (3e), i. e. , A mn , 

B m n* and C mn j van i s h except when n = 0 and m = v ^ 1. In other 
words, the scattered field consists of only two plane waves, propa- 
gating in two directions determined by the direction of incidence and 
the wavelength-to-surface period. These waves are both in the plane 
of incidence (i. e. , x-z plane) , and the scattering angle is determined 
by the equation sin 0 g = sin 0^ - X/L*. It is important to note that these 
waves exist physically only when | sin 0 S ] = | sin 0^ - \/ L | < 1 ; when 
\/L is such that this inequality is not satisfied, the propagation number, 
b(m, n) is pure imaginary, indicating an evanescent or non- radiating 
mode. For period much larger than wavelength, scattering takes place 
close to the specular direction. As \/L becomes larger, the two 
scattered plane waves diverge from the specular direction until they 
reach grazing ( | sin 0 g | = 1) , at which point the plane waves are cut 

off. Hence the scattering directions are determined not by the ampli- 
tude of the sinusoidal surface (so long as it satisfies the ''slightly rough 1 ' 
criterion) but by its period compared to wavelength. This behavior is 
also exhibited by the random rough surface considered below. 

Before discussing the scattering, let us discuss the parameters 
involved in the statistical description of the surface, ^(x, y) . Davenport 
and Root (1958) employ Rice T s earlier work (1944, 1945) to show that 
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a random but periodic function of period L can be expanded in a 
Fourier series such as Eq . (2), where the coefficients P(m, n) are 
random variables. Going further, they show that these coefficients 
become uncorrelated in the limit as L — * oo . Practically, it is 
necessary only that L >> l, i.e., the total surface length must be 
much greater than the roughness correlation length; this is exactly 
restriction (d) above. Now, recalling that the mean value of the sur- 
face is zero, these conditions become. 


( 6 a) <P(m, n)> = 0 


( 6 b) 


<P(m, n) P(u, v)> = 


f 0 for U) v = m, n 



W(p,q) for u, v = m, n 


where p = am = 2irm/ L and q = an = 2tt n/L. The function W(p, q) is 
physically the roughness spectral density of the surface, and p, q are 
radian wave numbers. Using these equations, the following relation- 
ships are obtained. 


(7a) 


<£ 2 (x, y)> = £ <P(m, n) P(u,v)> e iax(m " u) + ia ^ (n " v) 


m fi n,u, v 


- / <P(m, n) P(-m, -n)> — • - \ dm j dn W(p> ( l) 


-00 -00 


m , n 


1 r _ 2 

= 4 J j W(p, q) dpdq = h , 


- -00 


(7b) <Ux,y)t(x\y')> = £ <P(m. n) P(u, v )> + ’ aux '+ (any + iavy' 

m,n,u,v 

V \ rst iam(x-x') + ian(y-y ! ) 

= 7 <P(m, n) P(-m, -n)> e v w 7 

m,n 

^ 4 f ^ W(p, q) e P X ^ y dpdq = h 2 R(T x , T y ) 

* / 

-oo -o o 
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where t x = x-x’, and Ty = y-y 1 * The quantity h 2 , as mentioned earlier, 
is the surface mean square height, and R(t x , Ty) is the surface height 
correlation coefficient* Relationship (7b) merely states that the rough- 
ness height spectral density and the surface height correlation function 
are Fourier transforms. 

Now, let us determine the horizontal (h) and vertical (v) com- 
ponents of the plane wave scattered field in a given direction represented 
by the angles B s , cj> s ; first it is necessary to relate the wave numbers 
am and an of the exponentials of Eqs . (4b)-(4d) to the propagation con- 
stant of the plane waves in spherical coordinates . This relationship 
is obvious: 

(8) am = k Q sin 9 S cos 0 S , and an = k Q sin B s sin cj> s . 


Then the horizontal and vertical components of the scattered plane 
wave in the direction corresponding to m, n are 


(9) 

and 

(9b) 


tps / \ / A , „ x ia(mx + ny) + ib(m, n)z 

•E h (m, n) = (-^Vnn sm 9s + Bmn cos cj> s ) e 7 , 


E v (m, n) = (A mn cos 0 S 
- C mn sin 0 S ) 


cos <{> s + B mn cos 0 S sin <|> s 
g ia(mx + ny) + ib(m, n)z 


Let us redefine these relationships in terms of quantities 
( a hh anc * a vh) which are directly proportional to scattering matrix 
elements, namely 


(10a) 

and 

(10b) 


ia(mx + ny) + ib(m, n)z 


Eh(m, n) = 2k Q Ey a hh cos 0^ P(m-v, n) e 


q , a t->/ \ ia(mx + ny) + ib(m, n)z 

E® (m, n) = 2k 0 Ey a vh cos P(m-v, n) e y 


The quantities and a vh are determinable from Eqs. (9) and the 

expressions for -Amn* B mn » E mn - They will be given in simplified form 
subsequently, along with expressions for the other two elements a^v and 


a for vertical incident polarization. Note that the random nature of 
vv 

the rough surface is reflected explicitly in the surface coefficients 
above, P(m-v, n). The factors and a v h are nonrandom, known 

functions. At this point, the average scattering cross sections can be 
determined. First, let us find the increment in average field intensity 
in all of the plane waves between mi and m 2 , ni and n 2 . For arbitrary 
polarizations, y> this is given by 


m 2 n 2 


( 11 ) 


A < Ey (m, n) > = 4k 0 


E? 2 cos 2 6^ 



a ^6 1 

m ! TX\ 

<P(m-v,n) P*(m-v, n)> 


where use is made of the orthogonal properties of the coefficients, as 
stated in Eq. (6b). If the length of the surface, L, is large compared 
to wavelength, X. , there are many such plane waves between mi and 
m 2 (as well as between ni and n 2 ). The terms in the summation are 
then nearly constant for mi < m < m 2 (as well as ni < n < n 2 ) for m 2 
sufficiently near mi, and the summation can be rewritten 


( 12 ) 


A < | E® (m, n) j 2 > = 4k Q 2 I E 5 | 2 cos 2 Qi | tty ft 


<P(m-v, n) P#(m-v, n)> (m 2 -mi) (n 2 -ni) 


4k o h l E ?l 2cos20 i i ° v - 


W(am- aw an)AmAn, 


where Am = m 2 - mi and An = n 2 - ni • 

Using the above expression, the increment of intensity per unit 
of solid angle, S2 s , is expressed in terms of the Jacobian, AmAn/Afls • 


(13) 


A<|e^(b s , Ts )| 2 > , 4k 0 , 2 | E o|, cos2o . la YE J 2 


A 


L2 


AmA n 


_ v\i , , i aim 

W[ k Q (sin U s cos 4> s - sin 1^), k Q sin u s sin 4> S J — 
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Noting that Q s = sin b s dfr s dcf) s , the Jacobian is 


(14) 


AmAn 
A £2 s 



cos 8 S 


Hence the average field intensity at the receiver, or observation point, 
is the quantity in Eq. (13) multiplied by the solid angle subtended by the 
surface from the observation point. This is 


A^ s = 


L cos 0, 
R 2 


so that 


(15) 


,2 A < 

E v s l > = 


E° (0 S ,4> S ) 




L 2 cos 0 S 
R2 


The average scattering cross section per unit surface area at the 
observation point and for incident, scattered polarization states 5, y is 
given by 


or 


( y s><t>s) 


< 1 E§ 1 Z > 4ttR 2 

L2 fEg | 2 “ 


(16) o -°5 (0 S , <t> s ) - 4-rrko 4 cos 2 Bi cos 2 y s I a yD I 2 

W[k Q (sin 0 S cos (|> s - cos Oj), k Q sin 0 S sin 4> s J . 


Restriction (c) assumed that the surface roughness is isotropic. 
This means that the height correlation function R(t x , Ty) is a function 
only of the separation, r, between the surface points x, y and x\y', 

(i.e., r = J(x-x') 2 + (y-y'H = Jt x 2 + T y 2 ) . Hence, define 
( 1? a) p(r) = R( t x * T y ), 

applicable when the roughness is isotropic. This means that the rough- 
ness spectral density, W(p, q), is also symmetrical in p and q. Define 
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t = ^Jp 2 + q 2 and 

(17b) w(t) = W(p, q). 


Then for an isotropic surface the Fourier transformation between the 
roughness spectral density and the correlation function, as given in 
Eq. (7b) reduces to the following Hankel transform pair. 


(18a) 

and 

4h 2 p(r) 

= Ztt 

r p> 

\ tw(t)J Q (rt) dt, 
o 

(18b) 

w(t) = 

2h 2 , 

TT 

\ rp( r) J 0 (rt) dr . 

‘■'o 

Then Eq. 

. (16) may be rewritten 

(19a) 


*s> = 

— h 0 4 h 2 cos 2 3 i cos 2 0 S 1 'i v ,- | 2 • 1 

t r * ' 

where 




(19b) 

1 = W (D = 

r° J 

■ Z it \ r p (r) J c ( rt)dr , 
o 

(19c) 

t = k G > 

j sin 2 

0^-2 sin 0^ sin 0 S cos cj> s + sin 2 0s 


When the average incoherent scattering cross section per unit 
surface area is written as in Eq. (19a), the effect of the surface rough- 
ness is contained entirely in the factor h 2 I. The dependence upon 
incidence and scattering geometries, polarization states, and surface 
material properties is contained in the factors cos 2 cos 2 
As with all objects whose size is small compared to wavelength, the 
scattering cross section varies inversely as the fourth power of wave- 
length . 


Let us now interpret Eqs . (19) physically. They indicate that the 
average intensity of the scattered field in a given direction varies indirect 
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proportion to the surface roughness spectral strength at surface rough- 
ness frequency, t (radians / meter) . The quantity, t, in turn is a function 
of the incidence and scattering directions, as given in Eq. (19c). For 
example, for scattering near the specular direction, t ~ 0, the scatter- 
ing cross section is directly proportional to the roughness spectral 
strength at and near DC. On the other hand, the highest roughness 
spectral frequencies which can affect the average incoherent scatter- 
ed power are those at t = 2k 0 . This occurs where Qj = = *n7 2 and 

4> s = tt, i.e., for backscattering near grazing. Hence, the only possible 
surface roughness spectral frequencies which affect the average scatter- 
ed power in any arbitrary direction occur in the range 0 < t < 2k 0 ; any 
higher roughness spectral frequencies present in the surface can have 
no effect on the scattered power, at least to the first order perturbation 
theory used here and subject to the restrictions listed above* The 
physical interpretation of the scattering process is similar to that for 
the sinusoidal surface mentioned previously. 

Let us now examine two familar choices for the surface rough- 
ness correlation coefficient, p(r), and find the resulting spectral 
densities w(t) and the corresponding quantity I. 

(1) Gaussian Surface Height Correlation Coefficient: 
p(r) = exp ( - r 2 ij[ z ) . 

Here, i is referred to as the surface height correlation length. 

The quantity I for this case becomes 

(20) I a Tr J e 2 exp(-t 2 £ 2 /4) 

= tt£ 2 exp [ - (k Q 2 £ V 4 ) x (sin 2 0^ - 2sin 0^ sinO s cos cp s + sin 2 0 g )] . 

Surfaces having such a Gaussian surface height correlation coefficient 
are smoothly curving with derivatives at all points . The total mean 
square slope of the surface at any point with the above correlation 
coefficient is 


( 21 ) 



(2) Exponential Surface Height Correlation Coefficient: 
p(r) = exp (- | r | / 1 ) • 
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The quantity I in this case becomes: 


Ztt l 2 2 tt | 2 

( 22 ) 1 = = 

[1 +t 2 i 2 ]V 2 [1 +ko 2 | 2 (sin 2 0i-2 sin 0^sin 0 g + sin 2 0 S )J 3/2 

Surfaces with exponential correlation coefficients are jagged and 
have many vertical facets. The mean surface slopes and all higher sur- 
face derivatives are infinite, obviously because of these vertical facets. 
As an example, an urban area having buildings and houses comprising 
its scattering surface has a surface height correlation coefficient which 
behaves like the exponential, especially near the origin. 

It should be noted that for Zk Q i<l, the behavior of I for the two 
correlation coefficients is not significantly different, except for a 
factor of Z. Nor are the above two correlation coefficient models 
necessarily the most applicable in all cases. Peake (I960) found that 
a model which empirically fitted measured points for the correlation 
coefficient of an asphalt surface was 


p(r) = [ 1 + 20 r z f' Z 


where r is in cm. By the Hankel transformation between I and p(r), 
in this case I is an exponential function of t. 


C . Results for Vertical and Horizontal 
Polarization States 


We obtained above the average incoherent scattering cross section 
for a horizontally polarized incident plane wave and for the horizontally 
and vertically polarized scattered components. We left them in terms 
of scattering matrix elements, :< ^ and a v h» which are determined by 
Eqs . (9) and (10). The resulting cross sections and scattering matrix 
elements, a an d n vv> f°r a vertically polarized incident plane wave 
may be found immediately by noting that a vertically polarized plane 
wave (vertical or $ directed E-field) is identical to a horizontally 
directed H-field vector. Hence the quantities a hv an d ° vv niay be 
obtained from o v h and a hh by inte rchanging the roles of the relative 
constitutive constants, € r and p r , and employing care in observing 
the resulting signs of the E-field vectors from a knowledge of the H- 
field vectors. (The right subscript represents the polarization state 
of the incident wave and the left subscript represents the desired scat- 
tered wave . ) 
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1 . Bistatic scattering matrix elements 
for homogeneous surface 

(23 a) a hh = 

[ Mr cose i + > J«rH.r’ s ^ n2e l] [ fi r cos0 s + J « r p r -sin 2 e s ] 

€ r (fi r -l) J^rMr-sin 2 ^ - |i r (€ r -l) 

(23b) cl v h “ sin c() s . 

[fir cose i + N | € r» x r-sin 2 e i J [c r cosB s + j€ r p. r -sin 2 0 S ] 

(23c) ahv= sin <t> s • 

(c r - 1 )(€ r sin0jsin0 s - cosc(>s J^rM-r - sin 2 0^ s ^ n2 0 s ) + €r 2 (y*r~ 1 ) cos 4>s 

(23 d) ot vv = 


In the above expressions, * r and (i. r , the relative surface material 
permittivity and permeability, may each be either real or complex 
(indicating lossless or lossy material), and be either greater, equal to, 
or less than unity in magnitude. 

2 . Backscattering matrix elements for 
homogeneous surface 

The above elements are readily reduced to the case of back- 
scattering where <\> s = it and 0 S = . 


[.€ r cos0i + j€ r |j. r -sin 2 0i] [ « r cose s + ^€ r jJL r -sin 2 0 S J 




l)((i r sinBisin0 s -cos<t> s > « r |x r -sin 2 0^ J€ r ji r -sin‘ 


9 s) + Hr 2 


(€ r -l)cos<t> s 


(24 a) 

(24b) 

and 


^ _ (ji r - 1 ) [(H-r"l) sin 2 0^ + CrMr] ~ Mr 2 (*r“^ 

hh [ cos 0 { + Jc r f! r - Sin2 e t J 2 

a vh = O hv = 0 > 


(€ r -l) [(* r “M sin 2 + CrM-r] “ € r 2 (M-r - ^) 
[< r cos 0^ + J« r |Jir ~ sin 2 fc*i] 2 
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3 . Bistatic scattering matrix elements 
for perfectly conducting surface 


The scattering matrix elements for a perfectly conducting surface 
may be readily determined from Eqs . (23) by putting |j. r = 1 and € r — ► a>; 
this is evident from consideration of the basic boundary conditions for 
the two cases * 


(25a) 

a hh = - cos ’ 

(25b) 

sin <t> q 

«vh = - ’ 

COS Og 

(25c) 

sin <)> s 

P hv - + a 

cos 

(25d) 

sin B- sin - cos cj) s 

a vv cos Bj. cos Bg 


4 . Backscatte ring matrix elements for 
perfectly conducting surface 

Here again, we set cf> s “ tt and G s = B- \ n the above equations* 


(26a) ° hh “ ^ » 

(26b) o vh = a hv = ° , 

1 + sin 2 dj 

(26c) a vv = cos 2 0 . 


Actually, a perfectly conducting surface implies that the imaginary 
part of € r approaches infinity; the same result is obtained for either 
limit, however. 
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The average incoherent scattering cross section per unit surface 
area is now obtained using Eqs. (19) along with the proper scattering 
matrix elements, for the d ^ s ^ red combination of vertical and 

horizontal polarization states. These are given above. The factor I 
used in Eq. (19a) is given in terms of the surface roughness spectral 
density, as discussed previously and shown in Eq. (19b). 

Two outstanding facts are in evidence from the above equations 
for the important case of radar backs catte ring . For one thing, the 
cross polarized components are zero. It should be pointed out, how- 
ever, that these components, as indicated above, are zero only to the 
first order. The next order neglected terms in 0 *°^ and c are pro- 
portional to the mean square height and slope, k Q 2 h 2 and s 2 , with 
respect to c and cr° v . An expression for these terms is derived by 
Valenzuela (1967). The lowest order terms retained here vary with 
k Q 4 while the next order terms derived by Valenzuela are of order k Q 6 
in frequency. Thus our results indeed represent the solution in the low 
frequency limit. The second fact is that backscatter differs for the 
vertical and horizontal states. The vertical component is always 
greater than the horizontal for j €r I > |pr| » except at normal incidence 
where they are equal (i . e . , 0^ = 0). This fact has been observed experi- 
mentally, but previous theories based upon the less accurate tangent 
plane approximation have failed to show any such difference. 

(Davies, ibid.). 

Plots showing the average incoherent backscattering cross section 
per unit surface area are presented in Fig. 2 for perfectly conducting 
surfaces along with dielectric surfaces. The cross sections for both 
vertical and horizontal polarization states are given. These cross 
sections are normalized by dividing by ko 2 h 2 • The Gaussian correlation 
coefficient model given above was used with two relative correlation 
lengths (i.e., k c £ = 1.0 and k 0 £ = 5.0). The results for k 0 £ <1.0 are 
not significantly different in shape from those for k 0 £ = 1.0, but the 
magnitude of the cross section itself varies in direct proportion to k G 2 £ , 
as seen from Eqs. (19a), (20), and (22). Note that as correlation length 
(i.e., k 0 £ ) increases, more of the backscattered power is concentrated 
near the specular direction at normal incidence and less near grazing. 
This is physically reasonable, because surfaces with a given mean 
roughness height, h 2 , are smoother when the roughness height corre- 
lation length, £, is longer. This can be seen from Eq. (21) for the 
mean square slope of a rough surface. 
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Fig. 2. Average incoherent backs catte ring cross sections per unit 
area with linear polarization states from slightly rough 
surface model vs. incidence angle for Gaussian surface 
height correlation coefficients, normalized to ko 2 h . 

(1) k 0 £ = 1, vertical polarization states, (2) k 0 * = 
horizontal polarization states, (3) k 0 f =5.0, vertical 
polarization states, (4) k 0 i =5.0, horizontal polarization 
states . 
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D. Results for Circular Polarization States 


It is relatively easy to transform from the vertical and horizontal 
polarization states to the circular polarization states by using standard 
matrix methods. As mentioned previously, the elements a vv > a hh' 
a v h» an d a hv are directly proportional to scattering matrix elements, 
preserving both proper relative amplitudes and phase differences. 
Hence, a set of analogous elements applicable for circular polarization 
are determined and given below for use in Eq. (19a). 


1 . 


Bistatic scattering matrix elements 


* 


(27a) a LR = 


a hh + a w + U a hv - a vh) 


(27b) 


_ a hh “ a vv + U a hv "*■ a vh) 
a RR , 


(27c) a RL - 


a hh + a vv “ i( a hv “ a vh^ 


(27d) a LL 


a hh " a vv _ i( a hv j“ a yh) • 


The quantities on the right side of these equations are given in 
Eqs . (23) for a homogeneous surface material or by Eqs . (25) for a 
perfectly conducting surface. 


* The right subscript represents the polarization state of the incident 
wave and the left subscript represents that of the desired scattered 
wave . 
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2 . Backscattering matrix elements 


In the case of backscattering, it was shown in Eqs . (24) that 
a vh = a hv = 0- Hence the above elements reduce to the following: 


a hh + 

(28a) a LR - a RL = vv 


(28b) 


a LL = a RR 


a hh " a vv 
2 


In this case, the elements on the right side are given by Eqs* (24) for 
a homogeneous surface material and by Eq. (26) for a perfectly con- 
ducting surface. 

Physically, &RL and a,L R represent the ’’polarized” components 
in the scattered field; for backscattering from a smooth surface or 
symmetric object, the scattered circular state is always opposite 
sense from that of the incident circular state. ’'Depolarization 1 ’, then, 
appears in the elements QtLL and a RR> representing the scattered 
circular state of the same sense as the incident state. It is important 
to note that depolarization occurs for the circular states for back- 
scattering even though it does not occur for the linear states (i.e*, 
a hv = a vh - 0); this is due to the non- zero difference between 
and O'yv * 

Curves showing the normalized average backscattering cross 
section per unit surface area as a function of incidence angle are 
presented in Fig. 3. These represent perfectly conducting and dielectric 
surfaces for both the ’’polarized” and ’’depolarized” components. The 
Gaussian surface height correlation coefficient model is used with two 
values of the correlation length, k 0 £ = 1.0 and 5.0. It must be 
emphasized again that all of these curves were derived from a first- 
order perturbation theory; higher order terms in k Q h and s (relative 
rms roughness height and slope) have been neglected. 
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Fig. 3. Average incoherent backscattering cross sections per 

unit area with circular polarization states from slightly, 
rough surface model vs. incidence angle for Gaussian 
surface height correlation coefficient, normalized to 
k 0 2 h 2 . (1) kof = 1.0, opposite sense polarization states, 

(2) k G f = 1.0, same sense polarization states, 

(3) k c £ =5.0, opposite sense polarization states, 

(4) k o j0 =5.0, same sense polarization states. 
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E. Results for Arbitrary Linear 
Polarization States 


In many experimental situations, the transmitting and receiving 
antennas employ linear polarization states, but they may not necessarily 
be oriented along either the vertical or the horizontal polarization 
directions defined by the surface coordinate system. Define as the 
angle between the unit vector and the desired incident E-field com- 
ponent in the direction toward the unit vector . Define t] s in the same 

way with respect to 0 S and Then the matrix element a for 

bistatic scattering can be written s * 


( 29 ) 0L r\ s r] i = a vv cos r )i cos u s + a vh sin Ui c °s r\ s 

+ a^v cos 'Hi s ^ n + a hh s ^ n ^i s * n 'Hs 


1 . Backscattering with aligned and 
crossed linear polarized antenna 


Often the transmitter and receiver antennas are aligned so that 
“Hs = r li* For the receiver antenna polarization direction orthogonal to 
this, T) s = rji + tt/ 2 . Employing these angles and also the fact that 
a v h = a hv = 0 for backscattering, we arrive at the following two elements 

(30a) a riiTU = a vv c °s 2 "Hi + a hh sin 2 . 

and 

(30b) a ‘ r li+ir/2. T lx = '( a vv " a hh^ sin cos ^ ' 

These elements can be used in Eq. (19a) to determine the back- 
scattering cross section. The elements and cx vv to be used in the 
right side of the above are given in Eqs . (24) for a homogeneous 
surface material and Eq. (26) for a perfectly conducting surface. 
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2 . Angularly averaged backscattering 
with aligned and crossed linear 
polarized antennas 


In many cases aligned and crossed linear transmitting and 
receiving antennas may be employed, but the antenna angle, with 
respect to the surface vertical may be rotating and unknown at a given 
instant of time. This may arise in the case of satellite observations 
of a planetary surface. In other situations it can result from uncon- 
trollable Faraday rotation of polarization due to the ionosphere. An- 
other situation where an average quantity is desired is the case where 
a relatively short pulse of linear polarization sweeps past a planetary 
surface, illuminating an annular portion of the sphere. The angle, r^, 
between the linear direction and the local plane of incidence at a given 
point on the annulus varies continuously between 0 and 2tt in this case. 

In the above cases, the quantities desired are | a r|jr||| 2 and 
l Qr li+ 7 r/ 2 , r|^ | 2 » averaged by integrating over from 0 to 2tt and 
dividing by 2 tr . These quantities become 


(31a) < |a^ rj. I > - — [3|a.vvl + ^R e ( a vv a hh) + 3 | a ^ | ] > 

1 1 O 

and 


(31b) 


< l aT li+Tr/2, T li 


®vv - a hh I 
8 


The notation Re(x) denotes the real part of x. 


F . Comparison with Measured Results 

It is difficult to find measured data with which to compare the 
preceding theoretical results for a slightly rough surface for two reasons, 
(i) It is necessary that roughness height be everywhere small compared 
to wavelength, (ii) It is difficult experimentally to separate the average 
coherent scattered power from the incoherent; in most cases, no attempt 
is made to do this. For backscattering, the coherent contribution to 
the total power dominates near normal incidence. The results presented 
in this section apply only to the incoherent scattered power. 

Recently, average backscattering cross section measurements 
were reported by Wright (1966) from a fresh water surface at X-band. 
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Capillary waves were generated under controlled conditions such that 
the wave amplitudes were small compared to radar wavelength. Return 
for the both vertical and horizontal states were obtained, and are shown 
in Fig. 4 as a function of incidence angle. Along with these measured 
curves, theoretical curves are drawn, using the equations derived above 
and the parameters of the water surface estimated by Wright. These are 
c r - 55(1 + i.55), k Q h = .05, k Q i = 2. It should be noted that the measured 
results were presented not in terms of backscattering cross section per 
unit area, but backscattering cross section in square meters. Since the 
area of the illuminated surface was not given, the absolute value of cr° v 
and cannot be ascertained. Hence, the magnitudes shown in the 
figure are made to conform to the theoretical. The significant observa- 
tion is the agreement in functional dependence between the measured 
and theoretical curves. 

Another set of measurements on an asphalt surface at X- and 
K a *bands allow another comparison with the model. (See Fig. 5). 

Actual measurements of surface statisitics and dielectric constant 
were made in this case and used to compute the curves for the slightly 
rough surface model, (Cosgriff, I960). Again the functional dependence 
on frequency, angle, and polarization is in good agreement. Even the 
somewhat poorer agreement in absolute value is reasonable in view of 
the large number of surface parameters which had to be estimated 
from rather small samples. 


III. VERY ROUGH PLANAR SURFACE 
A. Introduction 


The remaining class of rough surfaces which can be treated 
analytically is the one whose members have a large scale roughness 
such that the local surface radii of curvature over nearly all surface 
regions are much greater than wavelength. In most cases of practical 
interest, the rms roughness height, h, will then also be larger than a 
wavelength* Then all of the scattered power from such a surface is 
incoherent, (i.e., the phase angle of the scattered field becomes 
uniformly distributed between 0 and 2tt, and hence the average value of 
the scattered field is zero; only the average scattered power is non- 
zero.) 

When this radius of curvature criterion is satisfied, the scatter- 
ing problem can be solved by several optics techniques. If these are 
used properly, they all lead to the same result. The value of examining 
more than one approach to the problem lies in the physical insight into 
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Fig . 4 . 


Measured average back- 
scattering cross section 
per unit area from slightly 
rough fresh water surface 
at X-band vs. incidence 
angle for vertical and hori- 
zontal polarization states . 
Dashed curves represent 
slightly rough surface model 
using parameters measured 
from fresh water surface . 




ANGLE FROM VE RT ICAL ( DEGREES ) ANGLE FROM V ERT ICAL ( DEGREES ) 


Fig. 5. 


Measured average backscatte ring cross section per 
unit area from slightly rough asphalt surface at X- 
and Ka-bands vs. incidence angle for linear polari- 
zation states . Dashed curve is slightly rough sur- 
face model using measured parameters p(r) = (1+Z0r 2 ) 
where r is in cm; k 0 £ == 0 . 085 at X-band and 0.30 at 
K a -band;k 0 i = ~ 0.48 at X-band and 1.7 at K a -band; 
€ r ^4.3+i0.1 at X-band and 2.5 + iO . 65 at K a -band . 


V 2 , 
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the scattering process which each provides. At present, three different 
optics approaches have appeared in the literature. Rather than retrace 
a derivation in detail here, a critique of these existing analyses will be 
made, along with the physical explanation of each. Finally, the common 
results of these analyses will be presented and curves for two models 
will be shown for comparison. 

All of the three optics approaches apply to the class of rough 
surfaces which satisfy (at least) the following restrictions: 

(a) The local surface radii of curvature at nearly every point 
on the surface are considerably larger than wavelength. 

(b) 3t,/3x, 3^/3 y <1.0, i.e., surface slopes are relatively 
small. Only when this is true can one neglect shadowing 
and multiple scattering over most of the range of incidence 
and scattering angles. 

(c) <(3£/3x) 2 > = <(3C/3y) 2> > i-e., the roughness is isotropic. 

The restriction is not necessary to obtain an answer, but is 
employed in this paper for simplicity. 

(d) L > > £ > > X, i.e., the dimensions of the illuminated area 
are much greater than the roughness correlation length, £, 
which in turn is much greater than wavelength. The second 
portion of the inequality is automatically fulfilled if (a) is 
true. The first part of the inequality precludes the possibility 
of looking at only a part of a hill or wave, but rather it de- 
mands that several hills or waves be included in the illuminated 
area . 

(e) k 0 2 h 2 >> 1, i.e., mean square roughness height, h 2 , is of 
the order of, or greater than, wavelength squared. 

(f) The surface height correlation coefficient, p(r), must be 
parabolic in behavior at the origin, i.e., p( r)^_^ 3 1 - r 2 / £ 2 , 
where r is the separation between two points on the surface 
and i is termed here the correlation length. 

The first two restrictions above permit use of the tangent plane 
(or physical optics) approximation. This approximation involves writing 
the scattered field at the surface as the product of the incident field and 
the proper reflection coefficient. Restriction (f) is imposed in order to 
avoid mathematical difficulties which lead to erroneous results when 
one employs a correlation coefficient whose behavior is linear at the 
origin, (Barrick, 1968). The second, third and fourth restrictions are 
identical to those given previously for the slightly rough surface. 
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Although the slightly rough surface analysis of the preceding part 
takes into account multiple scattering and shadowing, the normal optics 
techniques used for very rough surfaces ignore these multiple inter- 
action effects. Generally, they are higher-order effects so long as the 
surface slopes are not too precipitous and the scattering angles are not 
too close to grazing. No known correction factor for multiple scattering 
has been computed, but shadowing corrections have been formulated 
recently. Beckmann ) 1 965) first formulated the '’shadowing function.” 
This was later revised by Brockelman and Hagfors (1966) and Wagner 
(1967). The latter shows that the correction factor is less than 3 dB 
for angles of incidence and scattering more than 10° above grazing and 
an rms total surface slope, tan^s = tan" 1 (2h/f), of less than 27° . In 
view of the restrictions placed upon the class of surfaces considered 
here, shadowing corrections will not be introduced in this paper. 


B . Review and Interpretation 

1 . Physical optics approach 

The first approach to appear historically is also probably the most 
mathematically thorough and satisfying. By the same token, it probably 
sheds the least physical insight on the scatter mechanism because of the 
mathematical detail. This method employs a physical optics integral 
expression for the scattered field. The total fields at the surface 
appearing in the integrand are determined from the tangent plane approxi- 
mation mentioned above. The first comprehensive treatment of this 
problem for a perfectly conducting surface was made by Isakovich ( 1952) 
who employed a vector physical optics formulation. This work appears 
to have been overlooked by Western investigators, for Davies (1954) 
solves the same problem by the same techniques, but employs a less 
complete scalar formulation. Beginning about I960, a variety of 
Western and Soviet investigators began to extend these results. Among 
these are Hagfors (I960, 1964), Beckmann (1963), Hughes (1962), Fung 
(1964), Hayre (1961), Daniels (1961), Semenov (1965), and Stogryn (1967) 
Semenov appears to have been the first to have solved the problem for 
scattering from a surface of homogeneous material, € r ,p r , using a vector 
formulation, which accounts properly for polarization. The results of 
Semenov have been derived independently and simultaneously by Barrick, 
and the two are in agreement, after some algebraic simplification of 
Semenov’s results. Stogryn’s solution is the same as that of Semenov 
(although his notation is more complicated), and hence serves as a third 
check on Semenov's analysis. 

Semenov expresses the scattered vector field (for both vertical 
and horizontal incident polarization directions) in terms of a physical 
optics integral. For example, for an incident plane wave whose E-field 
is horizontally polarized, the integral, after application of the tangent 
plane approximation, can be written as follows: (See Fig. 1). 
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(32) 


E s = 


ik 0 e ik o R o 


4ir R, 


C r 


— \ - k s ] . r , 

F(^, r ) e oL 1 sJ dxdy , 

L/2 


where 


R, 


distance from origin to observation point, 


r = xx + y y + £,(x,y)z = distance from origin to point 

on rough surface, 

kp £ s - unit constant vectors pointing in direction of 
incidence and scattering, 

t» x , £y - local surface slopes in x and y directions at 
surface point £ (x, y), i . e . , 8£ / 3x and 3 £,/ 3y . 

The factor F(^ x , is a function of the local normal to the surface 
and of the local Fresnel reflection coefficients at each surface point. It 
can be expressed as follows: 


— >1 1 + ^x + r A A A — A A . Av— , 

(33a) F(i x ,i y )= — { n —\ . ly.tki^njGi + y.tkiCosVi + nJGJ, 


where 


(33b) G x = [cos Yi - cos Ys R j.(Yi)J L ( k i x n ) - k s *( k i x n ) k s] 

A A A 

+ [ 1 + RJ^)] kg X (ki + n cos Yi)> 

(33c) G u = [1 + R n (Vi )] [( k i + n cos Y\) * k s-( k i + " cos Yi) k s ] 
+ [cos Y{ ~ cos Y s R h (Yi)] [ k i cos y s - (lc s - ^ ) nj , 


A _ a A 

(33d) n = x ■ — = local unit normal vector to surface at 


Jl+^x + £ 2 point £ (x, y), 


y« , y s - local incidence and scattering angles in planes of incidence 
and scattering with respect to local surface normal, (i.e., 
cos y* = - ^ . n, cos y s = £ s .n), 
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R^Yj), R |( (-y^) - Fresnel reflection coefficients at local surface 
point t,(x, y) with incidence angle, Yi* 

As one can observe, the integrand consists of the factor F(C, x ,C,y), 
which varies with the local surface slopes from point to point on the 
surface, and an exponential phase factor containing the phase difference 
between the local surface point and the scattered field point. The first 
factor is to be removed from the integrand as a constant; the justification 
for this is the stationary phase (or specular point) argument. This 
argument maintains that for large k Q in the exponential factor, the only 
surface regions which contribute to the integral are those smoothly 
curving portions which are in a position to specularly reflect into the 
^esired scattering direction. The direction of the surface normal, 
ng p, and the local incidence angle, Yi = Ys = L » at these specular points 
are readily determined by bisecting the angle between the incidence 
propagation direction and the desired scattering propagation direction. 

The constant factor, F(£ x ,£ y ), is then evaluated using the slopes 

of the surface at these specular points. 

The remaining integral containing the exponential phase factor 
is not solved by the same specular point or stationary phase approach, 
however; the scattered field is squared to give scattered power or 
intensity, and then it is averaged. Averaging consists of multiplying 
the double integral by the joint probability density functions for the 
surface height at two different surface points and then integrating over 
these two random variables. The order of integration is interchanged, 
and the integrand of the former double integral is averaged first. This 
interchange is justified mathematically (See Middleton, I960) so long as 
restriction (f) is applied. With the aid of restriction (e), the remaining 
integral is solved in a straightforward manner; nearly all of the above 
mentioned references carry out this latter process correctly. 

The step of averaging under the integral sign, though mathe- 
matically correct, obscures the physical understanding of the scattering 
mechanism. It is also somewhat less than satisfying to justify the 
removal of a complicated factor in the integrand by the specular point 
or stationary phase argument, but then never to actually employ this 
stationary phase principle to solve the remaining integral. However, 
the two approaches to be discussed subsequently do justify the validity 
of this procedure by actually showing that the specular point theories 
give exactly the same result, at least for a Gaussian distributed surface 
height . 

No one has yet justified correctly a more exact method of evaluating 
this integral than the stationary phase approach. This is obviously due to 
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the extremely complicated form of the factor F(t >x ,C y) • At least two 
previous analyses have proposed alternate schemes for evaluating the 
integral. Barrick (1965) suggested averaging this factor independent 
of the exponential factor; hence, F(£x> £y)^would be removed as a con- 
stant evaluated at the mean slopes, i.e., F(0,0). Fung (1965) proposed 
expanding the factor in a series in £ x , £ y, and averaging each term 
separately. One can show that the resultant series cannot be truncated, 
since higher order terms are of the same order as the first terms. 

This is true even for vanishingly small values of surface slopes, £ x 
and Both of these approaches are inconsistent with the stationary 

phase principle, and therefore both can have no theoretical justification. 

Recently, Kivelson and Moszkowski (1965) have shown that for 
backscattering one can reduce the physical optics integral directly to 
the probability density function for the surface normal; this thus 
provides a valuable connection to the ray optics approach of the next 
section. Gaussian surface statistics are not necessary for their 
analysis, but the generalized form selected for their conditional sur- 
face height probability density function is still somewhat restrictive. 

The joint probability density function which they form is then not 
symmetric in the two height random variables, £ and £ \ as it must 
be, except for the special case of Gaussian statistics. 


2. Ray optics approach 

Recently, Muhleman (1964) brought an older - but better under- 
stood-theory to bear on the very rough surface problem. The "ray 
optics" approach (following his terminology) provides a more readily 
and simply explainable interpretation of the scattering process, even 
though mathematically it is less exact and elegant. The rough surface 
is initially approximated by a grid of small flat planar elements, all of 
which are connected to form the rough surface. Each element reflects 
incident power specularly, and the direction of reflection from a given 
facet is determined by the direction of its normal. Hence the amount 
of power scattered into a given direction is directly proportional to the 
number of facets whose normals are pointed in the proper direction 
(i.e., the direction bisecting the angle between the incidence and scat- 
tering directions) divided by the total number of facets. This quantity 
is the probability density function for the surface normal. It is 
expressible in terms of the probability density function for the surface 
slopes, since the surface slopes are perpendicular to the local normal. 

Muhleman dealt with the problem of bistatic scattering from a 
perfectly reflecting surface, but the problem of the homogeneous surface 
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material is just as easily solved. Recall that in the latter case, the 
components of the local E-field incident upon a planar facet in and 
perpendicular to the plane of incidence each have different Fresnel 
reflection coefficients, R u and R x . Referring to Fig. 1, we can 
determine the vertical and horizontal reflection coefficients for a 
properly oriented facet which reflects power incident from direction 
Q{ into direction B s ,(b s . Vertical and horizontal here refer to the 
polarization directions of incident and scattered fields with respect 
to the mean surface, i.e., the x-y surface, and not with respect to 
the local surface facet. Let us define R v h as the reflection coefficient 
between the horizontally polarized incident (right subscript) and the 
vertically polarized reflected (left subscript) E-field components with 
mean incidence and scattering angles ®s> $s • The other three co- 
efficients are defined in the same manner. Then assuming that scat- 
tering in these directions is produced by planar facets properly 
oriented to specularly reflect, these vertical and horizontal reflection 
coefficients are determined from straightforward geometrical con- 
siderations. They are: 


(34a) 


Rhh = 


sin 0-i sin B s sin 2 <t> s R M ( L ) + &2 a 3 R x ( l ) 
4 sin 2 i cos 2 i 


(34b) 


R v h = - si* 1 4>s 


a 2 sin 0 S R n (0 - a 3 sin 0^ R^(0 
4 sin 2 i cos 2 l 


(34 C ) 


R hv = " sin 


a 2 sin 0 S R x (0 - a 3 sin d\ R M («-) 
4 sin 2 L cos 2 i 


(34d) 


R 


sin 0-i sin 0 S sin 2 4> s R x ( i) - a 3 sin 0i R (| (i ) 


vv 


4 sin 2 t cos 2 i 


where 


(35a) 

and 

(35b) 


R,i ( T i) = 


Rx(^l) = 


f 

7 . 

€ T cos ^ - S 

€ r Hr “ sin *1 

e r cos t| + 

|«rM-r - sin 2 r\ 

|J. r COS T| - ^ 

€ r (j. r - sin 2 r| 

jx r cos 11 - 

i € rl JL r - sin 2 'H 
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and where l is the local angle of incidence at the specular reflecting 
facets, and is given by 

(36) cos 2 i = — ( 1 - sin sin B s cos c|> s + cos cos G s ) • 


The quantities az and a 3 involve the scattering angles and are given 
subsequently in Eqs . (54). As a check, the above reflection coefficients 
for backscattering reduce to the following: 


(37) Rhh* ^vv ~ R||(0), * R x (®)i ^hv* ^vh 


Since polarization differences were not considered by Muhleman, 
the above quantities squared must be inserted as factors in his results 
to obtain power scattering for the vertical and horizontal polarization 
directions as defined in Fig- 1. 

Hagfors (1966) showed that the probability density function for 
surface slopes appearing in Muhleman's result for scattered power 
can be easily related to the surface height probability density function 
when the surfaces are Gaussian. In doing so, he showed that for 
backscattering, one obtains the same result as from the physical 
optics analysis discussed above- It is readily shown that the two 
methods are equivalent for bistatic scattering also. [Note, it was 
pointed out by Hagfors that due to an oversight, one must divide 
Muhleman's result by the cosine of the angle between the local surface 
normal and the mean plane normal, i.e., by 


(38) 


a a 

cos (n, z) = 


cos 0^ + cos U s 
2 cos i 


• ] 


The ray optics approach sheds much needed insight into the 
scattering process from a rough surface. As a rigorous derivation, 
however, the method can be seriously questioned. Such an approach 
considers only scattered energy or power, and not the fields themselves 
hence phase relationships between fields from different regions are 
ignored from the outset. The ray optics method assumes specular 
reflection from a small planar facet in the form of a non-diverging 
tube or parallel rays. Later when one permits the planar facet to 
shrink to zero, the paradox between ray optics and modern geometrical 
optics theory becomes obvious. Diffraction effects and patterns, due to 
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phase interference, and divergence of rays due to curved phase fronts 
are ignored. The method, nonetheless, serves our purpose here be- 
cause, under the restrictions mentioned previously, we are dealing 
with incoherent scattered power. This means that it is sufficient to 
add power scattered from different portions of the surface, rather 
than the fields themselves . Hence ray optics is justified for this 
problem, and provides needed insight into the scattering process. 


3 . Geometrical optics or stationary 
phase approach 


The third approach is the application of the stationary phase 
principle to the Kirchoff integral for the complex scattered field. 

The result shows that scattering from a portion of a quadric curving 
surface does indeed radiate specularly. The scattering cross sections 
of such a curved surface is 


(39) cr = ir |Ri R 2 | , 


where Ri R 2 is the product of the two principal radii of curvature at 
the specular or stationary phase point. A very rough surface con- 
sists of many such specular points. Everyone has observed the 
dancing specular points on a rippling lake surface on a moonlight night. 
However, neither of the previous theories directly relate the total 
scattered power to the number of specular points and their radii of 
curvature, even though the physical optics method justifies the 
simplification of the integrand on the basis that this stationary phase 
principle is valid. It is therefore instructive to examine the scattering 
problem formulated on this basis and note that one arrives at the same 
result as the preceding two approaches. In addition, this third 
technique provides further insight by explicitly deriving the average 
number of specular points on a rough surface and their average 
Gaussian curvature, Ri R 2 • 

Kodis (1966) formulated the rough surface problem rigorously 
in this manner and showed by stationary phase that the scattering cross 
section for a rough surface can be expressed, as one intuitively expects, 
as follows: 

N 

/ ^ I II 1 ? i(4>i - 4 m) 

(40) ^ = V Z ( I R li R 2i I • I R ij R zj I ) e 

j =1 
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where the total number of specular points is N and the Gaussian cur- 
vature at the i-th points is RiiRzi- The quantity - 4>j contains the 
path length difference between the i-th and j-th specular points. 
Employing the restriction p, h > \ , one can show that for a very rough 
surface, the specular points are randoinly and widely separated in 
terms of wavelength, and therefore, and 4>j are uniformly distri- 
buted. Hence the average scattering cross section becomes 

N 

(41) < 0 -> = IT < > |RiiR 2 i I > • 

— / 

i-1 

Assuming now that one can assign an average Gaussian curvature to 
the specular points, one has 


(42) 


<a> = it <N> <! R] R z i> 


which one expects for incoherent scattering. The average scattering 
cross section per unit area, c r°. for incident and scattered polarization 
states - and y, can now be expressed 


(43) 0 - tt n < 1 R, Rz! > 

where 


(44) 


n 


<N> 

T7“ 


is the average number of specular points per unit area and Ry^ is the 
reflection coefficient for the surface at a specular point. For the 
vertical and horizontal polarization states (i.e., y> ” v, h), these 
coefficients are given in Eqs . (3Z). 

Bar rick (1968) used the above formulation and obtained expressions 
for n and <|RiR 2 (>. These results are identical to those of the previous 
two approaches. While confirming the general validity of and similarity 
among the various high frequency optics techniques applied rough sur- 
faces this third approach also provides valuable insight into the scat- 
tering process, which is complementary to that of the other two. In 
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particular, one finds that the average number of specular points is 
proportional to the probability density function for the surface slopes 
evaluated at the required slopes at the specular point; this quantity 
is greatest where the local normal at the specular points (or bisector 
between incidence and scattering directions) is near the mean normal, 
or vertical* This is reasonable from consideration of the moonlit lake 
example. The mean value of the radii of curvature product at the 
specular points, on the other hand, becomes greater as the local 
normal is farther removed from the mean normal* 


C. Results 

1 . General form 

Since the solutions for all three approaches are essentially the 
same for the Gaussian surface, we shall present one result here. 

It is expressed in a form which proceeds naturally from the physical 
optics approach and is obtained from the derivation of Semenov (1965). 
The average scattering cross section per unit surface area, for 
incident polarization state § and scattered state y, is written as follows: 


(45) 




J . 


The quantity (3^6 has exactly the same significance as 0.^5 for the 
slightly rough surface model: it is directly proportional to the scatter- 
ing matrix element relating the incident field of polarization state 5 
to the scattered field component of polarization state y. These quantities 
will be given below for the vertical and horizontal as well as circular 
states . 


The factor J is related to the roughness statistics in the following 
manner: 


(46) 




r J Q ( tr ) M££i(iu, -iu;r) dr. 


The quantity t is a function of the scattering angles and is given in 
Eq. (19c). The symbol u, also a function of these angles, is defined as 


( 47 ) u = k Q (cos 0 ^ + cos B s ) . 
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The function i(iu, iv;r) is the characteristic function of the surface 
height random variables £ and measured at surface points (x, y) 
and x'.y 1 ). (Formally, it is the double Fourier transform of the joint 
probability density function, P(C»£ * ; r); r is the horizontal distance 
between the surface points; 

(48) r = > (x-x') 2 + (y-y 1 ) 2 . 


Therefore, in order to obtain a closed form result for the scattering 
cross section, one must either know or assume a form for the surface 
height joint probability density function. (For the slightly rough sur- 
face, recall, one has to choose a form for the surface height correlation 
coefficient.) We shall select two common probability density function 
models for the sake of comparison and find the quantity J. 


a . Gaussian surface height 
joint probability density 

Using the standard Gaussian joint probability density function 


(49) 


p g (£’C) 



[?> 2 £'p(r)+ £' 2 ] 

2h 2 [ 1-p* (r)] 


along with restriction (f) on the form of the form of the correlation 
coefficient p(r) near the origin and restriction (e) on the size of the 
mean square roughness height, h 2 , compared to wavelength, we can 
readily solve Eq. (46) for J. 

J_ t 2 

(50) J = _±_ e " s 2 ^ , 

s 2 u 2 

where 



is the total mean square slope of the rough surface. 
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b . Exponential surface height 
joint probability density 


For the sake of comparison with the preceding model, we select 
an exponential probability density model of the following form: 


(51) 


Pe(C. V) 


3 

2 it h 2 [ 1 -p 2 (r)]2 


* ¥- 2 iV P ( r )+^ 2 * 
J h [ 1- p 2 (r)] 

V- ^ 


1 

2 


Equation (46) for J is solved using the same restrictions as with the 
preceding model. 

n/ 6 t_ 

(51) J = -1?— e ' S U 

S 2 U 2 


The optics techniques show that for very rough surfaces, the 
only parameter reflecting the degree of roughness in the scattering 
cross section is s, the rms roughness slope. The rms roughness 
height, h, does not appear explicitly, and as such is not easily 
determined from observations of average scattered power from such 
a surface . 


2. Polarization dependence 

The elements (3^6 now be given for several polarization 

states . 


a. Vertical and horizontal states 

For the vertical and horizontal polarization states, the scattering 

matrix elements (3w> Phh> Pvh* and (3hv ma Y be found upon simplifi- 
cation of Semenov's results, or from either of the other two approaches 

to the problem. They are: 


(53a) 


vv 


sin sin 9 g s in 2 4> s R x (i) + a 2 a s R (| (t) 
a i a* 


(53b) 


Phv = ' sin 


a 2 sin 0 S R ± (t) - a, sin 0^ R n (0 

a l a 4 
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(53c) 


Pvh = ' sin <t>s 


a 2 sin 0 S R x (l) - a* sin 0* R„(l) 


a i a 4 


(53d) 


Phh = 


sin Qj sin B s s jn 2 <fr s R^i) + a 2 a s RJi) 
a i a 4 


where R h ( rj), R x ( "n) and i are defined previously in Eqs . (35) and (36) 
and ai, a 2 > a 5 , a* are given below. 


(54a) 

(54b) 

(54c) 

(54d) 


ai = 1 + sin sin 0 S cos <\>q - cos 0^ cos 0 S , 

a 2 = cos 0^ sin 0 S + sin 0^ cos 0 S cos 4> s , 

a 5 = sin 0{ cos 0 S + cos 0^ sin 0 S cos 4> s , 

a 4 = cos 0^ + cos 0 S 


From the above expressions and the two probability density 
models previously given, one can substitute the P's and J into Eq . (45) 
to obtain the average scattering cross section per unit area. 


b . Circular states and arbitrary 
polarized linear states 


The P's given above for the vertical and horizontal states, being 
directly proportional to the scattering matrix elements, can be sub- 
stituted into Eqs. (27) of the preceding part to give the elements PlR> 
PRR’ PRL» and PLL ^ or c ^ rcu l ar polarization states. Equations (29) 
and (30) of that section can be used to give Pr^r) s for linear polarization 
states arbitrarily directed. 


3 . Results for backscattering 

Backscattering is the most important and commonly occurring 
situation for the radar engineer. As such, the above results will be 
specialized to this case, and curves for the two probability density 
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models as a function of incidence angle will be presented. The expres- 
sions for the (3's given in Eqs. (53) are allowed to approach the back- 
scattering limit (i . e . , Q s — ' * ®i» 4>s it ) to yield 


(55a) Pvv~ Phh ~~ “ sec 6 i R„ (°) = sec 9 i R ± (0) , 

(55b) Pvh = Phv 0 • 


The circular polarization backscattering matrix elements found 
from substituting these expressions into Eqs. (27) are 


(55c) P LR = p RL “ P v v = Phh sec °i R„ (°) * 

(55d) PLL = PRR 0 * 

The backscattering matrix element between arbitrary linear 
states as defined and given in Eq. (28) becomes 


(55 e ) pT^Tig = Pvv cos ( T ls - “Hi) = " sec 0^ R h (0) cos(rj s - r} { ) . 


The above equations demonstrate the claim of Hagfors and other 
investigators who note that optics techniques applied to very rough 
surfaces models predict no depolarization for backscattering. The 
second and third optics models discussed previously clearly show that 
backscattering can come only from areas or facets on the rough surface 
which are oriented normal to the incidence direction. This also 
accounts for the appearance of the Fresnel reflection coefficients for 
normal incidence in the above equations . 

Let us now give the results for average backscattering cross 
sections per unit area for the two probability density models proposed 
above, remembering that no depolarization is predicted. 
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a. 


Gaussian surface height probability 
function model 


( 56 ) 


vv 


0 -° = <r? r = <r° = (T? p = <T° 

hh rt, LR r)-r)| 
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sec 4 ©. 
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b . Exponential surface height 
probability function model 


(57) 
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vv 
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hh 


= cr 
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o 4 

3 sec b- 
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tan 6- 


The first of the above equations has been given previously by 
Hagfors (1964) while the second is presented here for comparison. 
Curves illustrating the dependence of 


| r (! (°)| 2 

for the above two models upon angle of incidence, for various values 

of rms roughness slope, s, are presented in Fig. 6. 


IV. COMPOSITE SURFACES 
A. Explanation 

Natural forces seldom create a "very rough" surface satisfying 
the restrictions of the preceding part. Rough surfaces arising in nature 
having roughness whose height is large compared to wavelength almost 
always also possess smaller scale roughnesses, such that the surface 
is not locally smooth. As examples, consider an area of mountainous 
terrain. The mountains and valleys alone comprise a roughness large 





Angle From Verticol , degrees Angle From Vertical , degrees 

Fig. 6. Average incoherent backscattering cross section per 
unit area for very rough surface statistical model 
vs. incidence angle for any set of direct polarization 
states for two surface height JPDF (Joint Probability 
Density Function) models . Curves are normalized to 
the Fresnel reflection coefficients for normal incidence. 


compared to wavelength for frequencies of VHF or higher. If no other 
smaller roughness were present, the analysis of the preceding section 
might be applied. However, the trees, rocks, and grassy fields also 
present constitute a roughness whose scale compared to wavelength 
may be small. Hence, the tangent plane restriction and optics approxi- 
mation of the preceding part cannot be used without closer examination. 
The type of surface described above, characterized by both large and 
small scale roughnesses, will be called a "composite" surface. 

Let us look again at the chief differences between the average 
incoherent backscattered power from the slightly rough surface and 
the very rough surface of the preceding parts. The slightly rough 
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surface, as shown in Fig. 2, produces a small backscattered field near 
normal incidence where 0^ = 0, but this field may still be significant 
near grazing. As discussed previously the scattered power near graz- 
ing is directly related to the amount of surface roughness spectral 
strength at higher spatial frequencies; in particular, at and near 
frequency 2k 0 radians/meter. On the other hand, backscattering 
from very rough surfaces, as seen in Fig. 6, is much larger near 
normal incidence (0i = 0) but falls off rapidly away from normal. 

The suggestion to be put forth here is that both the large scale 
and the small scale roughness contribute significantly to the scattering 
from natural, composite surfaces. Near the specular direction (i.e., 

0s = ^i* 4* s = 0) scattering of the optics type, predicted by the very 
rough surface theories, predominates. At scattering angles consider- 
ably removed from this direction, however the smaller but persistent 
scattered power is due to the slight roughness present on top of the 
larger roughnesses. 

In fact, to a first approximation, one may merely add the 
average incoherent scattering cross sections from the very rough 
surface model to that of the slightly rough surface model to obtain 
the total composite rough surface scattering cross section. An 
heuristic and physically intuitive proof of this statement will be given 
here, but for a more exact-albeit mathematically involved- proof, one 
can consult two recent Soviet articles of Semenov (1966) and Fuks (1966). 

Consider first the very rough surface model of the preceding 
part. The surface to which it applies is smooth over regions of the 
order of a wavelength. The scattered power is incoherent and con- 
centrated near the specular direction because the surface slopes are 
generally small and the probability is great that the local normal of 
any surface point is within a few degrees of the vertical. Now assume 
that a slight roughness "grows" on top of this very rough but locally 
smooth surface. This slight roughness does not greatly diminish the 
large specular scattering of the very rough surface, just as it did not 
diminish greatly the coherent scatter of the smooth plane over which 
it was imposed in Part 2. However, in directions considerably away 
from the specular (e.g., near grazing for backscatter), where the 
very rough surface model predicts very small scattered power, the 
slight roughness will produce an observable addition to the scattered 
power. Since the "very rough" portion of the surface has small sur- 
face slopes, the slight roughness can be approximated to a first order 
as being distributed over a smooth plane in computing the local angle 
of incidence. As seen in Figs. 2 and 3, small changes (say 15® or 
less) in the local angle of incidence, 0^> (due to local variations in the 
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underlying very rough surface slope) do not produce a marked change 
in backscattering cross section, so long as the spectral strength of 
the superimposed slight roughness is fairly constant out to approxi- 
mately 2 k Q radians/meter. Hence the two average incoherent scat- 
tering cross sections may be simply added together, assuming the 
same incidence and scattering angles (d^, 6 S , 4>s) f° r both. 

Actually, two averaging processes must take place. The slight 
roughness height variable and the large roughness height variable are 
averaged separately. This necessitates the additional restriction that 
the two different scales of roughness are statistically independent. 

Such is generally true when the two roughness scales are formed by 
different natural processes. For example, mountains and valleys 
are usually formed by an entirely different process than is the 
vegetation cover; hence the assumption of statistical independence 
between the distribution of these natural roughnesses, although not 
perfect, generally seems reasonable. Since the slight roughness and 
large roughness are statistically independent and each alone produces 
strictly an incoherent scattered power, we can perform the averaging 
of each separately and merely add the powers from each. 

The above explanation is meant to be intuitive rather than 
rigorous. It is meant to provide physical insight rather than mathe- 
matical confusion. Mathematical proof is satisfying, however, and 
can be found in the above mentioned Soviet references . Generally 
the necessary restrictions are those given in the preceding two parts, 
along with the one mentioned in the preceding paragraphs. Neglected 
in this explanation is any possible roughness scale which is of inter- 
mediate size such that it does not fall entirely either in the very 
rough or the slightly rough category. Since most natural surfaces 
have a roughness component whose scale is of this order, the model 
here must be expected to give results which deviate somewhat from 
measured values. This deviation is expected for scattering geometries 
where the scattering cross section for the slight roughness alone is 
approximately equal to that for the large scale roughness. In general 
the predicted result will be too small, because any intermediate scale 
roughness will tend to add another scattered field component which 
will be significant in this region. Since no known theory is valid for 
intermediate scale roughness, this component cannot at present be 
taken into account quantitatively- 

Western writers have attempted to treat composite rough sur- 
faces. Beckmann (1965) and Hayre and Kaufman (1965) both treat this 
problem by a physical optics model. Their analysis does not truly 
apply to a slight roughness scale since it is based on an optics formulation. 
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This is evident from the results, which are functionally the same as 
for the very rough surface with a single large scale of roughness. The 
only difference is that the mean square slope, s 2 , appearing in their 
results for the composite surface is written in terms of the sum of the 
mean square slopes of each of the roughness scales present. Their 
results are hence valid for all roughness scales down to those which 
can no longer be treated by their optics formulation (i.e., which no 
longer satisfy the restrictions of Part 3). They cannot apply to sur- 
faces which include a true slight roughness component, however. 

Radar experimentalists have in the past attempted to break the 
measured direct polarized backscatter cross section vs. incidence 
angle into two regions. The region near normal was called the 
specular region, and it was generally agreed to be predicted by the 
optics theories and results of Part 3. The "tail" of the echo, or that 
region near grazing was termed the "diffuse" region, and there has 
been no general agreement to date on a satisfactory explanation for 
its existence. Various investigators have attempted to fit empirical 
curves to this component; among the postulated are cos 0-p cos 3 / 2 d{, 
and still others. Hagfors (1967) has suggested a model consisting of 
randomly oriented dipoles as accounting for this component. This 
is one of a large class of geometrical models with distributions of 
individual scatterers which lie outside the class of continuous sur- 
faces we have been considering. (See Peake, 1967). 

We wish to suggest here that this "diffuse" component is pro- 
duced by the presence of slight and intermediate roughness scales 
and is predicted to a first order by the models of Part 2. We offer 
predicted curves based on the composite model presented here and 
compared them with measured results for backscatter. 

The explanation for the depolarized backscattered component 
from composite surfaces is more complex, however. The optics 
theories predict no return. To the first order, the slightly rough 
model studied here also indicates no depolarized return for the verti- 
cal and horizontal states; depolarization for the circular states does 
appear, but is absent at normal incidence. The higher order terms, 
for example the second term as calculated by Valenzuela, do show 
depolarization, even for the vertical and horizontal states. As the 
roughness scale increases with respect to wavelength from slight to 
intermediate, these higher order terms containing depolarization 
effects become non- negligible . Since the perturbation theory breaks 
down as the smallness parameters approach unity, the contribution 
of the intermediate scale of roughness cannot be calculated. It is clear 
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from the trend of the higher order terms, however, that they will con- 
tribute significantly to both the polarized and depolarized components. 

It is our conclusions, therefore, that depolarization for backscatter is 
produced by two mechanisms: (a) the slight and intermediate roughness 
scales just discussed, and (b) multiple reflections from suitably 
oriented surface elements of the large scale roughness for surfaces 
where the reflection coefficients are polarization dependent. 

To illustrate the predictions of the composite surface model, 

Figs. 7 provide curves for horizontal and vertical polarization states 
for a variety of dielectric constants and model parameters. They are 
made by simply adding the backscattering cross sections for the slightly 
rough model to those of the very rough model. More curves of this type 
can be found in Radar Cross Section Handbook (1968). Figure 8 shows 
the measured cross section of the sea surface on a calm day made 
for the vertical and horizontal states; these are also taken from the 
above reference. The similarity is convincing. Note that none of 
these curves represents a depolarized component; the received 
polarization state is in each case the same as the transmitted. 

For the circular polarization states, one can construct the 
same type of curves, which results are sort of a superposition of the 
curves of Fig. 3 with those of Fig. 6. They are not shown here but 
may be found in the preceding reference. It should be noted, however, 
that the curve representing the depolarized component is not actually 
a composite, because the very rough surface contribution to back- 
scattering in this case zero: only the slight roughness produces de- 
polarization. One set of curves for circular polarization is compared 
to a measured set from the lunar surface at \ Q = 68 cm in Fig. 9. 
Measured results were taken from the work of Evans and Pettingill 
(1963), and were normalized with respect to the measured cross section 
of the entire moon at this frequency. The set of model curves was 
selected on the basis of best fit. From the model for this best fit, the 
following parameters for the lunar surface are indicated: c r ~ dielectric 
constant =2. 9, tan^S] = large scale roughness slope = 12°, h 2 = rms 
slight roughness height ~ 2 cm, tan _1 s 2 = slight roughness slope = 27° • 
The comparison of the measured and predicted depolarized components 
shows lack of agreement near normal incidence. This is expected for 
the following reasons: (i) Hagfors (1967) reports that the ability to iso- 
late one measured circular state from the other was probably not 
greater than 18 dB . Hence near normal where the polarized component 
is very strong, most of the “measured" depolarized component may 
actually be the polarized component, (ii) The intermediate scale rough- 
ness, which cannot be accounted for by the theory, most likely produces 
depolarization. As mentioned previously, this effect is likely to appear 
further away from grazing, i.e., near normal incidence. 
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Fig. 7b. Exponential large scale roughness 
height probability density model. 







Fig. 8. Measured average back- 
scattering cross section 
per unit area from re- 
latively calm sea surface 
atX-bandvs. incidence 
angle for vertical and 
horizontal polarization 
states . 



V. ROUGH SPHERICAL SURFACE - 

BACKSCATTERING CROSS SECTION 

A. Introduction 


The past sections have derived the average incoherent scattering 
cross section per unit area for a rough planar surface, i.e., a surface 
whose mean value is a plane. In this section we will extend these 
results to a rough spherical surface. Such a situation is representa- 
tive, for example, of a planetary surface, or of many spherical 
satellites of the Echo I and II variety whose surfaces are somewhat 
rough. In defining the average backscattering cross section of a rough 
spherical surface, we understand that the entire sphere is immersed 
simultaneously in the illuminating field. This may not be the case, for 
example, in many of the planetary radar studies where a short pulse 
sweeps past the planet, illuminating only an annular area at a time. 
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Fig. 9. Measured average back- 
scattering cross section 
per unit area from lunar 
surface at UHF (X. = 68 cm) 
vs. incidence angle with 
circular polarization states . 
Dashed curves represent 
composite rough surface 
model, presented for com- 
parison. Upper curves 
represent opposite sense 
polarization states, lower 
curves - same sense. 




The results of the preceding sections can be used to obtain the 
incoherent backscattering cross sections for both very rough and slightly 
rough spheres. Furthermore, for the slightly rough sphere, there will 
be a significant coherent component due predominantly to reflection from 
the front cap, and this will also be estimated. 

It is assumed that the radius of the sphere Ar, is much larger 
than both the wavelength, X, and roughness height correlation length, 
i . It is further assumed that the sphere is large enough and its 
material lossy enough that waves entering the sphere do not further 
contribute to the scattering. In addition, all of the restrictions of the 
preceding sections apply. 
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B • Coherent Backscattering Cross Section 


It was mentioned in Part I that the coherent scatter from any 
slightly rough surface is highly dependent upon the shape of the illumi- 
nated surface in the absence of roughness. The presence of roughness 
is taken into account by modifying the Fresnel reflection coefficients at 
the surface, as in Eq. (1). The reflected complex fields in this case, 
being coherent and having a non- zero average value, can be treated 
just as they would be from a smooth surface; stationary phase analysis 
predicts that nearly the entire coherent backscattered field in this case 
comes from the front specular cap. Hence the average coherent back- 
scatter cross section from a slightly rough spherical surface becomes 

2 2 

(58) 0 - = it I R(0) | 2 e~ k ° h 
where 

N/ir 

(59) R(0) = — ^ _ 

N/€ r + N/H*r 



is the Fresnel reflection coefficient for a smooth planar surface at 
normal incidence* The effect of surface roughness near the specular 
cap is accounted for by the exponential factor involving' the mean square 
roughness height, h 2 . The above result, having come from the use of 
the slightly rough surface reflection coefficients of Eq. (1) is valid for 
surfaces having a Gaussian distributed roughness height. One sees 
that as rms roughness height increases to where k Q h > 1, the coherent 
cross section decreases rapidly. 

The polarization states of the scattered coherent field above is 
the same as it would be in the absence of roughness, i.e., there is no 
depolarization predicted. 


C . Average Incoherent Backscattering Cross 
Section - Slightly Rough Surface 


The average incoherent backscattering cross sections for a 
slightly rough spherical surface are computed from the preceding 
results for a rough planar surface. Use is made of the property of 
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incoherent power which permits one to merely add the power backseat - 
tered from one portion of the surface to that from another portion, 
rather than adding the complex scattered fields themselves . Recall 
that cr ^ is the average incoherent backscattered power per unit mean 

surface area (where 0 S = 6; and 4>s = tt in the equations for cr g). Then 
the backscattering cross section for a mean spherical area element 
dA = s in 0i d0i d<j> is the product of dA and cr . The total average 
backscattering cross section is then determined by integrating over the 
entire illuminated hemisphere, i.e., 


(60) 


V5 


2 

= A R 


Ztt 


f 


n / 2 


y5 


sin 0^ d0i 


where the integration over 4> has been performed since the integrand is 
not a function of cf> . The spherical coordinate system has been chosen 
here so that the polar axis coincides with the direction of propagation, 
making the polar angle 6^ also the angle of incidence upon the spherical 
surface . 

For the aligned and crossed incident and received linear polari- 
zation states, the expressions which must be substituted into the above 
equation for cr yf. are those given in Eqs . (19) and (31). The integration 
over <j> corresponds to the integration over which was already per- 
formed in deriving Eqs . (31)* The integration over 0^ demanded in the 
above equation can only be performed numerically. The resulting cross 
sections normalized as 


cr 

TT Ajf k Q 2 h 2 

are plotted in Fig. 10 vs k 0 ^ f° r various values of surface dielectric 
constant. This is done for both the aligned and crossed linear states 
(representing the polarized and depolarized return) and for the Gaussian 
and exponential surface height correlation coefficient models . Note that 
as k 0 l becomes larger the depolarized component becomes smaller. 

This is expected, since larger k 0 i corresponds to a smoother surface. 
Note that the terminology vertical and horizontal states have no signifi- 
cance when referred to a sphere. 

For the circular polarization states, the appropriate curves are 
found in Radar Cross Section Handbook (1968). They are nearly identical 
with those for the linear states in the cases of the polarized components j 


53 



Perfectly conducting surface 


IO.I3kft! — 


0.599tk|ll, 


€ r * 5.0 


€, * 2.0 


5.85 X lO 1 3 


I.ISxlo' 


€ r*'S 

^O.I205l4i 2 

'0.01478 


0.002226k*! 2 


1.25x10* 


Perfectly conducting surface 


2.133(4** 


€,• 5.0 


0.0265ll4!*- 


€ r « 2.0 


-0.001642k*!* 


-000003(75(4** 


0.0001397 - 

: ik 0 *) 4 j 

\ J 


Fig. 10a. Average incoherent backscattering cross section for 
slightly rough spherical surface of mean radians A r 
as a function of k D f- Cross section is normalized to 
it Ar k Q "h . Linear polarization states, upper curves 
represent same sense received as transmitted, lower 
curves - perpendicular sense received, (a) Gaussian 
height correlation coefficient model. 
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and plotted as a function of the arctangent of the rms roughness slope, 
s (i.e., s = (2h/l)). This function is sometimes referred to as the 
roughness gain of the sphere. The two curves represent the two 
probability density function models discussed in Part III. Noteworthy, 
however, is the fact that as roughness increases, the backscattering 
cross sections become larger than that for a smooth sphere. 

It is seen that as the surface becomes smooth (i.e., s — ^ 0 ) , the 
backscattering cross sections do indeed approach those for a smooth 
sphere i.e., cr = tt A 2 ^ |r( 0)| 2 . The figures are of questionable 
validity for tan -1 s > 45° because of the restrictions under which 
they were derived. For typical surface slopes much beyond 45°, 
shadowing and multiple scattering are no longer negligible. In the 
roughness slope region less than 45° the increase in scattering cross 
section from the smooth sphere limit is not appreciable. Hence an 
important conclusion for planetary surface backscatter can be drawn: 

Large scale surface roughness for all practical purposes 
can be neglected and the total average backscattering cross 
section can be taken to be that for a smooth sphere, 
cr = tt A 2 ^ | R(0) | 2 

Consequently, knowing the approximate radius of the body, Ar, an 
estimate of the surface dielectric constant can be made from a 
measurement of o' / tt Ar ; for the moon for the range 5 cm < k <100 cm 
one obtains approximately |r( 0)| 2 = 0.07. Such a value of | R( 0 ) | 2 
implies an effective dielectric constant of € r = 2. 9* 


E. Average Incoherent Backscattering Cross 
Section - Composite Surface 


As mentioned previously, most natural surfaces consists of 
roughness scales both larger and smaller than wavelength. As shown 
in Part IV, one can add the incoherent backscattering cross sections of 
the preceding sections (presented in Figs. 10 and 1 1 ) to obtain the 
effective total backscattering cross section of the composite surface. 
However, the backscattering cross section for the slight roughness, 
being proportional to k 0 2 h 2 , is small, since this parameter is re- 
stricted to small values. Hence only the large scale roughness contributes 
significantly to the backscattered power for the polarized received com- 
ponent. On the other hand, since only the slight roughness scale is 
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responsible for the depolarized component, according to these models, 
it alone may be the best available estimate for this component. 


VI. SUMMARY 

We have attempted to set forth what we believe to be the soundest 
analyses and solutions for rough surface scattering that are presently 
available. In addition, we have enumerated the restrictions under which 
the results are valid. Finally, we have interpreted physically the scat- 
tering mechanism behind the mathematics. 

A slightly rough surface produces strong coherent reflection 
near the specular direction and a weaker incoherent scatter in other 
directions . An estimate of the coherent reflection can be computed 
just as it is for a smooth surface, making use of the reflection coef- 
ficients of Eq. (1); the resulting scattering pattern is strongly de- 
pendent upon the shape of the illuminated area and exhibits a lobe 
structure. The incoherent scatter was dealt with in more detail by 
a perturbation technique rather than the less valid tangent plane approxi- 
mation. It was shown that incoherent scattered power from a slightly 
rough surface of either homogeneous material or a perfectly reflecting 
interface is directly proportional to the roughness spectral densities. 

The highest roughness spectral components which can affect the process 
are those near Zk Q radians/meters- These are responsible for back- 
scatter near grazing incidence. Lower roughness spectral frequencies 
produce scattering closer to the specular direction. For the important 
case of backscatter, significant differences exist between the cross 
sections near grazing for the vertical and horizontal states; vertical 
transmit and receive antennas produce considerably more return than 
horizontal, a fact which is confirmed experimentally. The solutions, 
including only the lowest order perturbation terms, exhibit no de- 
polarization for the vertical and horizontal states. From the solutions 
for vertical and horizontal polarizations, results for the circular states 
and arbitrarily oriented linear transmit and receive states are given. 

The only presently satisfactory methods for treating very rough 
surfaces are optics techniques; basic to all of these is the tangent 
plane restriction. A physical interpretation of the scattering process 
can be based upon any of the well known optics principles. Three 
techniques (i.e., physical optics, ray optics, and the specular point or 
stationary phase principle) all give identically the same result. Scatter- 
ing in a given direction takes place as a result of surface facets so oriented 
that they specularly reflect. As would be expected from such a theory, 
backsc attered power is not depolarized. Two probability density function 
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models are chosen for the surface height; for both, the curves of back- 
scattered power show strong return near normal incidence (the specular 
direction) which falls off rapidly as grazing is approached. The only 
surface height correlation coefficient which can be employed with these 
techniques is one which is parabolic in behavior for small surface 
separations . 

Natural surfaces which are very rough in height compared to 
wavelength almost always have smaller scale roughnesses superimposed; 
these are called composite surfaces. As one might expect, the slightly 
rough and very rough surface scattering theories together account in 
part for the return. In the case of backscatter near vertical, for 
instance, the return is dominated by the very rough surface predicted 
scatter. Near grazing, however, the smaller roughness scales usually 
account for the return. Curves for backscattering cross section are 
shown for these composite models at various set of parameters for the 
vertical and horizontal polarization states. These are compared with 
measured data from the lunar and sea surface. 

Backscatter from roughened spherical surfaces, which has 
application to planetary probing and passive satellite communitations, 
is treated in Part V. The coherent and incoherent average backscat- 
tering cross sections are obtained, the latter in graphical form resulting 
from a numerical integration. Very rough spherical surfaces are also 
treated, and it is shown that the roughness in this case has little effect 
on the backscatter relative to that of a smooth sphere; in reality, the 
roughness slightly increases the cross section. 
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